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ABSTRACT 

A three dimensional model of thermospheric dynamics is developed in terms of 
the eigenfunctions of the atmospheric system. Formulae for the external heat 
inputs like solar XUV- radiation and corpuscular heating during geomagnetic 
storms are derived in terms of these eigenfunctions. Those series contain tidal 
components (depending on local time) and planetary components (depending on 
seasonal time) which are the generators of tidal and planetary neutral atmospheric 
waves. The relative importance of corpuscular heating when compared with 
XUV- heating is estimated. Approximate analytic solutions for the generation and 
propagation of the atmospheric waves within the dissipative thermosphere, ex- 
cited by the solar heht input, are given. It is shown that the eigenfunctions which 
are the Hough-functions within the nondissipative lower atmosphere change into 
the spherical surface functions within the dissipative thermosphere. Moreover, 
the density amplitudes of the wave modes decrease proportional to l/n^ where n 
is the zonal wave domain number of the spherical harmonics. Therefore, only 
the wave modes with low wave domain numbers n are significant at thermo- 
spheric heights. 
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A THREE DIMENSIONAL MODEL OF THERMOSPHERIC DYNAMICS; 


I: HEAT INPUT AND EIGENFUNCTIONS 
I. INTRODUCTION 

The dynamic behavior of the neutral component of the upper atmosphere above 
about 100 km altitude has been studied extensively in the last decade by means 
of direct in situ measurements of the neutral density (e.g.; Taeusch et al., 

1968; Newton, 1970), the neutral wind (e.g.: Vasseur, 1969) or by satellite drag 
observations (e.g. Champion, 1970; Jacchia, 1970). These observations must be 
compared and confronted with variations in the ionized component of the thermo- 
sphere as inferred from ionospheric drifts (Harnischmacher and Rawer, 1968; 
Sprenger and Schmindler, 1967), meteor trail observations (Roper and Elford, 
1965; Greenhow and Lovell, 1960), and geomagnetic variations (Matushita, 1967), 
as well as from measurements of electron density (Rawer and Suchy, 1967) ion 
composition (von Zahn, 1970) and ion temperature (Mahajan, 1969); Waldteufel 
and McClure, 1969). The available data of the neutral component of the thermo- 
sphere reveal the following general picture of thermospheric dynamics: 

(a) The diurnal component strongly dominates the daily density variation at 
thermospheric heights above 200 km with relative amplitudes up to 0.5. The 

relativ'e amplitude and phase remain roughly constant with altitude above about 

•'ll 

300 km altitude. The semidiurnal component above 200 km is weak though not 
negligibly small when compared with the diurnal component; 
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(b) The meridional structure of the diurnal density variation during equinox 
follows roughly a simple sin 9 -■ law, where 6 is the co-latitude*, 

(c) The annual variation is relatively weak. Below 400 km altitude the density 
bulge follows the zenith angle (Jacchia and Slowey, 1967). On the other hand, 
there exists a relatively strong semiannual variation within the whole thermo- 
sphere which can be observed down to the stratosphere (e.g.: Cook, 1970; 

Cole, 1968); 

(d) The average neutral density as well as the diurnal density amplitude 
strongly depend on solar activity indicating the dominant influence of the solar 
XUV heat input on thermospheric dynamics; 

(e) During geomagnetic storms the density increases significantly suggesting 
additional heating of the thermosphere by solar corpuscular radiation and 
Joule heating. 

Some of these observations are consistent with ionospheric data. E.G., the dy- 
namic behavior of the ionospheric F-layer can partly be interpreted by means 
of neutral winds deduced from the neutral pressure field (Kohl and King, 1967). 
The wind field deduced from the geomagnetic Sq current at 115 km altitude shows 
a predominant diurnal component (Kato, 1956; Stening, 1969; Tarpley, 1970). 
However, some ionospheric observations are in apparent disagreement with 
measurements of the neutral component. E.G., meteor trail measurements and 
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ionospheric drift measurements at higher latitudes reveal the predominance of 
the semidiurnal wind component up to F-layer heights (Greenhow and Lovell, 
1960; Harnischmacher and Rawer, 19G8; Sprenger and Schmindler, 1967). The 
neutral temperature inferred from Thomson- scattering measurements gives 
maximum values at 16.00 local time which seems to be in disagreement with the 
density maximum at 14.00 local time as determined from satellite and rocket 
bom measurements (Mahajan, 1968). 

The neutral and ionized components of the thermosphere form a coupled dynamic 
system which reacts to the gravitational tidal forces of sun and moon, to the 
solar XUV- radiation, to solar corpuscular radiation and indirectly to the total 
solar radiation via dynamic coupling from the lower atmosphere. 

V 

\ theory of thermospheric dynamics should be able to interpret the above 
mentioned observations in a selfconsistent manner and to explain the discrep- 
ancies between the observations of the neutral and the ionized component. Such 
theory should lead to a clear separation between the various energy inputs and 
their influence on thermospheric dynamics depending on day, season and solar 
cycle. Because of the complexity of the problem we are far away from a final 
theory of such kind. 

Harris and Priester (1962) were the first to attack this problem. They calcu- 
lated the thermodynamic reaction of a vertical air column to the solar EUV heat 
input. Considering heat conduction as the most important energy transport. 
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process at thermospheric heights they were able to reproduce qualitatively the 
observed diurnal behavior of the neutr’al thermospheric density. However, 
quantitatively their results revealed a discrepancy in amplitude and phase be- 
tween theory and experiment which they resolved by introducing an unknown 
second heat source. 

Lagos and Mahoney (1967) extended this model to other latitudes and showed that 
the EUV heat input alone is not sufficient to create the observed large thermo- 
spheric temperatures at high latitudes. They concluded that horizontal energy 
transport or an other heating mechanism must be important for the heat balance 
of the thermosphere. Dickinson et al, (1968) calculated a two dimensional model 
along a zonal strip. They included horizontal longitudinal winds and showed 
that in their model the time of the diurnal density maximum shifted toward the 
early afternoon as the result of adiabatic heating due to vertical motions. 


A two dimensional equatorial model has been calculated by Volland and Mayr 
(1970). Here the radiation condition of characteristic waves within the thermo- 
sphere has been applied as botmdary condition. This allows the separation of the 
contribution from the various heat sources within and outside the thermosphere. 

It has been shown that within the lower thermosphere a tidal diurnal wave gener- 
ated below 120 km height dominates the dynamic features while above 200 km the 
thermospheric dynamics are excited primarily by the EUV-heat input within the 
thermosphere. The second heat source of Harris and Priester could be elimi- 
nated by the introduction of horizontal longitudinal winds. These winds are 
restrained by collisions with the ionospheric plasma. Heat advection due to these 
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winds shifts the phase of the density amplitude from 17 hours LT to between 14 
and 15 hours LT (depending on the numerical value of the collison frequency) 
consistent with the observations, and it reduces the density amplitude to the 
observed values. 

The two dimensional model of Volland and Mayr (1970) has been adopted to 
reproduce quantitatively the dependence of the density on solar activity (Volland, 
1969a). Likewise, amplitude and pliase of the geomagnetic activity effect, of the 
27-day rotation effect, and the semiannual effect have been calculated and show 
satisfactory agreement with the observations. 

A first approach to the next step - a threedimensional thermospheric model - 
has been made by Lindzen and Blake (Lindzen, 1970; Lindzen and Blake, 1970). 
They used the concept of equivalent plane gravity waves to study the generation 
and propagation of some of the more impoitant diurnal and semidiurnal tides 
within the thermosphere. 

The model of Lindzen and Blake includes the assumption that the latitudinal 
structure and the equivalent depth of the various tidal modes do not change with 
height within a dissipative atmosphere. However, it is well known (e.g. Siebert, 
1961) that the eigenfunctions of atmospheric dynamics on the rotating earth which 
are the Hough functions wi thin the lower atmosphei’e, become the spherical 
functions if one can neglect the Coriolis force. The neglect of the Coriolis force 
is indeed a reasonable approximation at F- layer heights and above where ion- 
neutral drag and the viscosity forces by far exceed the Coriolis force. 
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In this paper which consists of three parts we shall study various tidal and plane- 
tary wave modes, their latitudinal dependence, height structure and their genera- 
tion and propagation within thermospheric altitudes. In this first part we shall 
outline an analytical approach to describe these various modes with the aim of 
providing insight into some of their general characteristics, For that purpose a 
number of approximations are necessary. 

The main approximations adopted are the use of perturbation theory and the re- 
striction to gravity waves. It is well known that at thermospheric heights two 
other wave types — heat conduction waves and viscosity waves — can exist. They 
influence mainly the amplitudes of temperature and horizontal winds. However, 
at least the influence of heat conduction waves on thermospheric dynamics has 
already been studied in a two dimensional model (Volland and Mayr, 1970) and 
this indicates that in our simplified model the temperature amplitude is subject 
to errors which are not too serious and which should be tolerated considering 
other uncertainties involved and considering the simplifications gained. 

In this part, we shall furthermore develop the solar XUV and corpuscular heat 
inputs as well as Jacchia’s (1964) exospheric temperature distribution into series 
of spherical harmonics. This will show that only few components in this series 
with low wave domain numbers are required for a sufficiently accurate repre- 
sentation of the observed temperature distribution. It suggests that the spherical 
functions approximate rather well the eigenfunctions of the thermospheric dynamic 
system. Finally, we shall determine in this paper analytic solutions of the generation 


and propagation of the various eigenfunctions within an Isothermal lower 
atmosphere and within an isothermal thermosphere. 

The analytic solutions developed in this firfjt part shall be used in the second and 
third part of this paper to calculate explicitly the wave propagation of several im- 
portant tidal and planetary wave modes which are generated by the corresponding 
components of the solar heat input. These waves include the fundamental sym- 
metric diurnal and semi-diurnal modes, the antisymmetric diurnal mode and the 
annual and semiannual modes. In the second and third part of this paper we shall 
study moreover in detail the altitude variations in the latitudinal structure for these 
modes. Especially, it will be shown how the eigenfunctions which are in the Hough 
functions within the lower atmosphere transfer into the spherical functions at 
thermospheric heights. Furthermore a detailed analysis of the wind systems of 
the various wave modes will be given. 

A more sophisticated numerical study of the three dimensional thermospheric dy- 
namics including heat conduction waves and taking into account realistic tempera- 
ture profiles within the thermosphere will be given in an additional paper (Volland 
and Mayr, 1971a (referred to as paper I)). 

2. THE MODEL AND ITS RESTRICTIONS 
a. Application of perturbation theory 

In order to find tractable analytic solu tions for the three dimensional spherical 
model of thermospheric dynamics we Imve to restrict outselves to a model which 
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is as simple as possible. The most important approximation in our approach 
is the use of perturbation theory. That allows us to treat the various eigen- 
functions of the atmosphere as if they were decoupled from each other and to 
separate the thermospheric disturbances related to the various energy resources 
of different time scale. The usual method of numerical analysis of thermospheric 
data in fact implicitly adopts the same assumption when one resolves the various 
perturbation periods like diurnal, semidiurnal, annual periods ect. by means of 

some kind of spectral analysis. If coupling between the individual eigenfunctions 
would be significant, one had to expect a significant contribution from higher har- 
monics of the basic periods as well as mixed frequencies. This has not been 
observed. 

We can estimate the errors involv3d in such a perturbation treatment. The 
horizontal and vertical winds at thermospheric heights are of the order of 

!uj < 100 m/sec 

jw| 1 m/sec 

(Volland and Mayr, 1970). The total time derivative in the equations of mo- 
mentum and of energy conservation is for the diurnal tidal waves 

^ ^ + (v- grad) j fi + j k u + j k w 

Cl t o t ^ * 
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with 


k -^l.S X 10“'* km"^ (horizontal wave number) 


1 •• 2 ** 1 

k^ 10 km (vertical wave number) 

2 H„ 


= 7.27 X 10 ^ sec"^ (angular frequency of the earth's rotation) 


Hq is the pressure scale height. Thus, the nonlinear terms are of the order 


K i^i wi 


n 


n 


< 0 . 2 . 


They still allow a perturbation treatment if one accepts errors of the order of 
20%. If one considers the uncertainties in the experimental data as well as in 
our knowledge of the energy sources and the physical coefficients like 
heat conductivity, viscosity and ion- neutral collisions at thermospheric heights, 
the error due to the perturbation approximation seems tolerable. Moreover, it 
will be shown in the following sections that the thermosphere behaves like a low 
pass filter for tidal and planetary waves. Thus it suppresses waves with great 
wave domain numbers and therefore prevents the generation of higher harmonies 
due to nonlinear coupling. 

For the perturbated thermosphere the barometric height formula is assumed to 
be valid. That implies the neglect of the vertical inertial force which, when com- 
pared with the earth's acceleration force g leads to the inequality 

|w] « — 10® m/sec 

* ' r\ 
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valid for all tidal and planetary waves. Moreover, we shall neglect the horizontal 
viscosity forces and the horizontal heat fluxes which are small compared with 
the equivalent horizontal and vertical components respectively below 400 km: 


Tj 92 u/;» x2 


K rf2 T/'3 x2 


k2 

X 

r; ci 2 u/;3 2;2 


K ^2 T/Vi z2 


k2 

7 . 


(u = horizontal wind; T = temperature; r; = coefficient of viscosity; K = coefficient 
of heat conduction; is the earth's radius). Finally, we shall neglect the 
metric factors in the equations of mass and energy due to the spherical co- 
ordinates (r is the distance from the earth's center and z = r - Rq ): 


because of 


i- — ('r2. .'t - .. JL 

j -2 3 r 3 r fj z 



— 3x 10"2 


Perturbation theory presupposes the knowledge of the mean physical parameters 
of the steady state atmosphere. We shall adopt as such steady state atmosphere 
the Jacchia model (Jacchia, 1964) and consider the vertical structure of this model 
as the quiet thermosphere during a particular solar activity. All variations in 
time and space due to temporal and spatial variations of the heat sources are con- 
sidered as perturbations superimposed upon that quiet thermosphere. 
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We shall also neglect the temporal variations in the neutral composition. In the 
case of tidal waves such neglect leads to an underestimation of the temperature 
amplitude which however is not too serious (Volland, 1969a), In the case of some 
types of planetary waves the wind induced variations in the neutral composition may 
give rise to substantial errors at heights where atomaric oxygen predominates. 

This aspect of the thermosphere dynamics will be investigated in separate studies, 
b. Restriction to gravity waves 

Within the lower atmosphere, tidal waves are of the type of internal gravity 
waves. Within the thermosphere, dissipation effects due to heat conduction and 
molecular viscosity give rise to the generation of three new wave modes: heat 
conduction waves, and two kinds of viscosity waves (Volland, 1969b). These wave 
types become increasingly important with increasing height and contribute 
especially to the temperature and horizontal wind amplitudes of the tidal vari- 
ations. In this approach here we shall neglect those waves and restrict ourselves 
only to the gravity waves. However we shall introduce the effects of heat con- 
duction and viscosity in the propagation of gravity waves by replacing the second 
derivatives with height of temperature and horizontal winds by height dependent 
coefficients in the equations of horizontal momentum and energy conservation. 

That procedure in fact is equivalent to an outfiltering of the three additional 
wave types. 

In paper I we shall take into account the heat conduction waves in the numerical 
calculations. There, we shall discuss the implications due to the neglect of the 


additional wave types. 
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It will turn out that our approach in fact gives results which become only modi- 
fied by a more exact treatment. Temperature and horizontal wind amplitudes 
are most affected by our approach. However it is possible to simulate the in- 
fluence of heat conduction waves and viscosity waves by the introduction of a 
slightly modified vertical profile of the heat input and by an effective ion-neutral 
collision number respectively. 


The dissipation coefficients are determined in the following manner: From the 
work of Jacchia (1964) we know that at a given time the vertical temperature 
profile at thermospheric heights above 120 km can be rather well approximated 
by the formula 


T - T - CT ~ T 

^aj V *(jo ^ 120 ' 




( 1 ) 


Here T„ is the exospheric temperature, Tjjq is the temperature at z = 120 
km altitude, and H is a scale factor which is of the order of 


H “ 40 km 

during moderate solar activity. Assuming vertical profiles of the wave 
amplitudes in the form of equation (1), we obtain for the heat conduction term 
in the energy equation 


B z 



2Z,,c 


y pQ ni 
" 1) To 


( 2 ) 
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with 




2.5 (3T-T„) e 

4 7 T2 H2 


-( z- z 1 2U 1 H 


Here is p^ = the mean density; the mean temperature, the mean pres- 
sure, and y = Cp/c^-'-l.S the ratio between the specific heats, r; = A vT the 

coefficient of molecular viscosity; K = 2.5 0^77 the coefficient of molecular heat 
conductivity (Chapman and Cowling, 1952). It will be shown later that the wind 

profiles have approximately the same vertical structure as that expressed in ( 1 ). 

Therefore, we find for the viscosity force in the equation of horizontal momentum. 

with 

^ _ ■>! (U„ - Uj 2 q) (3 u - uj H 

4 fi u2 h 2 

where H is the scale factor of the wind profile. Provided H , both 
dissipation factors Z.,(, and are of the same order of magnitude. They 
tend to become constant with height because of 


T - T„, u and 00 exp (- z/H). 


An estimate using data from the Jacchia model gives 


^vi ! 


' 0.05 at z < 1 50 km 


0. 1 at z > 200 km 
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From the more exact numerical study it is inferred that the dissipation factors 
Z|,e and Zyis are indeed of the order of 0.1 to 0.5 above 200 km altitude. 

We conclude that molecular viscosity and heat conduction can be neglected 
when compared with the Coriolis parameter 1/2 below 150 km. However, these 
terms can play a significant role in thermospheric dynamics at heights above 
200 km. In this connection it is worth mentioning that the ion-neutral drag term 

Z , r JL (4) 

(v‘ is the collision number of ion-neutral collisions) is of the order of three at 
F2-layer heights (Dalgarno, 1964). Therefore, in this height range and for tidal 
and planetary waves, ion drag is generally more important than molecular 
viscosity (Geisler, 1967; Kohl and King, 1967). 

3. THE DISTRIBUTION OF THE ENERGY INPUT 

The thermosphere reacts to solar heating or to gravity forces and behaves like 
an oscillator system in which forced oscillations are generated. The external 
energy generator may have any possible spatial and temporal distribution. If 
our assumption about the validity of perturbation theory is correct, we should 
be able to develop the energy generator in terms of the eigenfunctions of the 
thermospheric dymmic system and to study separately the response of the 
individual thermospheric wave modes to the corresponding energy modes. 
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It is generally accepted that the influence of the gravitational forces of moon 
and sun on the atmospheric tides are small when compared with the solar heat 
input due to solar radiation (Chapman and Lindzen, 1970). Therefore, we shall 
neglect that excitation mechanism in the following. However, electric fields may 
be a significant source of wave generation. These fields originating either from 
the Sq region at E layer heights or from the magnetosphere are mapped up or 
dovm along the geomagnetic field lines and give rise to ion drift across the geo- 
magnetic field at F2 layer heights. This ion drift in turn leads to wind and pres- 
sure perturbations in the neutral component via ion-neutral drag (Hines, 1965a). 

The electric field of the Sq region is generated by the neutral tidal wind at E 
layer heights. It thus acts like a coupling link between lower and upper ther- 
mosphere. It couples a fraction of the kinetic energy of the wind field at E layer 
heights to the winds at F2 layer heights. Since the wind at E layer heights itself 
is generated by solar heating within and below that region, such coupling can be 
considered as a vertical redistribution of the total solar heat input within the 
atmosphere. In a similar sense the electric fields from the magnetosphere are 
responsible for the transformation of a fraction of solar corpuscular energy 
stored within the magnetosphere into wave energy of the neutral thermosphere. 
They generate electric currents within the ionosphere which in turn heat the 
thermosphere via Joule heating. They also force the ions to drift across the 
geomagnetic field lines and induce wind and pressure variations at F layer 
heights. 
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The solar heat input into the thermosphere via photo-dissociation, ion recombi- 
nation, and electron cooling is itself a very complicated process and far from 
being well understood. A very simple approach v*hich is probably not too far 
from the reality is to assume that the effective heat input per volume is pro- 
portional to the mean pressure at thermospheric heights. This assumption 
may even include the virtual heat input due to the electric fields. In this paper 
we adopt that picture with a slight modification. We take a height dependence of 
the heat input Q like 


Q 

Po 


“ g (z) 


'2h 


( 5 ) 


with a constant value of e near, but not necessary equal, to one. 100 km 
is a reference height above ground. For convenience we shall consider in this 
paper an isothermal thermosphere, = const. In the numerical study of 
paper I we shall use the value e = 1 and a realistic temperature profile of the 
thermosphere. 


It remains to determine the temporal and spatial distribution of the heat source. 
As basic period we select the time of one year with the angular frequency of 


n - — = 2 X 10"^ sec"^ (6) 

one year 

Higher periods, e.g. the eleven-year-cycle, are slow enough to be treated quasi- 
stationary. 
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We first consider the XUV input of the sun. In order to obtain an approximate 
distribution of the XUV input we assume that its effective heat input is propor- 
tional to the zenith angle y of the sun 


(1 + e) g (z 


<^xuv 


2 Ho - .-0^ 


cos Y 


( 7 ) 


for lyi 1 7f/2 , and zero otherwise. Here is the total height, space and 

time integrated averaged heat input (in erg/sec) above the reference height 

Zq = r^ “ Rq where r^ is the distance from the earth's center and is the earth's 

radius. Though the real distribution may deviate more or less from equation (7) 

it should nevertheless give the right order of magnitude in particular since we 

shall deal here primarily with the relative magnitudes of the various frequency 

components. 


We develop the function (7) into a series of spherical functions. Approximating 
the time of sun rise and sun set by 


- arc cos (tg 6 ctg 0) v/2 - tg fc ctg 6' 


(8) 


and taking 


cos 1 


sin ''' S 


where * = 5^^ cos t is the solar declination with a maximum value of = 0.4, 
t is the universal time (t = 0 at the begin of the year) and d is the polar distance, 
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we obtain the following series: 


^XUV ' 6 (z) 


p« - 2 s„ CCS n, t p« bj cos 2 n. t} p« 


(2PI - ^ «o ‘ ny^’> - 


s 2 T f , 


( 9 ) 


^xuv “ ) Qxuv ^ ^0 erg/cm3 sec) is the average heat input per 

volume at the height z. P™ (0) are the spherical functions in Schmidt's normali- 
zation depending on polar distance 0 , r = nt+X. is the local time and \ is the 
geographic longitude. 

Though the development equation (8) breaks down near the poles, the coefficients 
of equation (9) approximate the energy distribution of (7) with an accuracy of 
about ±10% which has been tested by an exact numerical calculation. In equation 
(9) all terms with second or higher powers of 5^ have been neglected (except 
those that produce the semiannual component). These neglected terms have 
magnitudes of the order of 0.1 or smaller. 

Turning to the solar corpuscular heating we consider only the mean heat input 
average over the individual substorms. The heat input due to precipitating fast 
electrons and electric fields from the magnetosphere is confined to the auroral 
ovals, and most of the heating occurs on the night time hemisphere where the 
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component of the polar electrojet flows (Akasofu, 1968). It is well known that the 
geomagnetic disturbances related to that heat input do not show a significant 
annual variation. However, the semiannual component with a maximum during 
the equinox is pronounced (Chapman and Bartels, 1951). The semiannual vari- 
ation of the u, - measure of Bartels can be represented by 

Uj fi (1 - 0. 1 cos 2 t) (10) 


where u depends on solar activify, Uj is a measure of the disturbed geo- 
magnetic horizontal component. Thus, the magnetic energy of the disturbed geo- 
magnetic field is proportional to uJ . Assuming proporlionality between the 
magnetic energy and the dissipated heat input, we arrive at a heat source aver- 
aged over storm time and space of 


^corp 2 fig t) (in ei g/sec). (11) 

Assuming that this heat input is deposited within an infinitely thin strip at ±65° 
latitude on the night time hemisphere, we develop it into spherical harmonics 
(Volland and Mayr, 1971b). 


<co rp 


- qo 

^corp 


g (z) (1 - 0.2 cos 2fi t) {po + 3,66 P" 


2. 22 PO 


(12) 


+ (0.81 P} + 3.58 Pj +. . .) cos T + . . .} 
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with 


(1 + s') 0 

Cfi _ ^ ' ^corp 

'*corp *’ ” 


877H 


the average heat input per volume at = r - Here, we used for convenience 
the same height dependence for Q as for . The higher order coefficients 
in equation (12) may, of course, change significantly if we choose any other lati- 
tude and they would decrease if we would use a strip of finite broadness. In any 
case, these coefficients already indicate that this series converges very slowly, 
and that we need a great number of terms to approximate sufficiently well our 
assumed heat input distribution. However, as will be seen in the next section, 
the atmospheric dynamic system filters out wave modes with great wave domain 
numbers (n, m) so that in fact only the energy modes of lower degree in equation 
(12) can excite significant density variations at thermospheric heights. 


It should be mentioned here, that this is true only for the storm time averaged 
corpuscular ^eat input. Any individual polar substorm may generate short 
periodic gravity waves which travel from the auroral zones into the lower 
latitudes (Testud, 1970; Chimonas and Hines, 1970). However, these waves are 
a local phenomenon and outside the subject of our paper. 

A third possible energy input into the thermosphere may have its origin within 
the lower atmosphere. Such energy input is due to wave energy dissipation of 
waves generated within the lower atmosphere and penetrating into the dissipative 
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thermosphere. These waves may include the whole spectrum of short periodic 
internal gravity waves (Hines, 1965b), the fundamental diurnal tidal wave (Vol- 
land, 1969a) or semidiurnal tidal waves (Lindzen and Blake, 1970b). Wave dis- 
sipation is a nonlinear process and therefore can not be treated in our pertur- 
bation theory. This wave dissipation affects essentially the zero component P° 
of the heat input and gives rise to an increase of the mean temperature. We 
therefore have to introduce an additional heat source 

0dis^^Qdisa-0-05cos2n, t)P0 (13) 

which contains a semiannual component. The numerical value of this semiannual 
component is anticipated from the following results. During moderate solar 
activity it is 


0 ° > 0 ° 

^dis - '^xuv 

indicating that a significant fraction of the average exospheric temperature is 
caused by that heat input due to wave dissipation (see Page 3 of this paper). 

Before we shall discuss the implications of equations (9), (12) and (13) on thermo- 
spher c dynamics, we want to develop also the exospheric temperature 
distribution of the Jacchia model (Jacchia and Slowey, 1967) in terms of 
spherical harmonics. Jacchia’s model is supposed to represent well the density 
variations above about 250 km. The result is (derived for a solar activity factor 
of F = 125 and the Jacchia parameters n’ = m' = 2.5) 



- T° [PO - 0,048 cos t P» « 0.007 P» 

+ (0.12 P} ~ 0.008 cos t Pp cos ('- - 15.0'’) 

+ 0.035 P2 cos [2 (t - 13.2*’)] - ] 


( 14 ) 


TO 

0 


- 0.04 PO 




cos 2 0 t 

a 


with TO = 1000®K. Here, the first bracket on the right hand side in equation (14) 
is an excellent approximation of the analytical function of Jacchia's exospheric 
temperature distribution while the second bracket in equation (14) takes account 
of the semiannual variation of which is eliminated in the Jacchia model. We 
related the observed semiannual temperature amplitude to the (O.O)-component 
in equation (14) and indicate our ignorance about any latitudinal dependence by 
terms with arbitrary coefficients c^ , which should exist due to the correspond- 
ing terms in equations (12) and (13) though according to Cook (1970) the semi- 
annual variation does not show a significant latitudinal dependence. 

One aim of this paper is to relate the individual terms in equation (14) to the 
corresponding energy components in equations (9), (12) and (13). That is 
not possible without some ambiguity because we observe in fact the 
thermospheric response of the combined energy sources. We shall return to 
these questions in some detail in the two following parts of this paper. However, 
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some general conclusions can already be drawn from a comparison of equations 
(9), (12), (13), and (14). 


First, the XUV heat source contains a rather strong zonal (2.0) component which, 
als already Lagos and Mahoney (1967) pointed out, should give rise to a very low 
mean exospheric temperature at higher latitudes. The temperature distribution 
of equation (14), on the contrary, shows a rather weak mean latitudinal variation 
with the temperature gradient directed towards the poles. Jaccia’s first model 
(with his parameters n' = 2.5 and m’ = 1.5) distinguished by an elongated pressure 
bulge has even a mean temperature gradient directed toward the equator equivalent 
to a positive (2.0)-component in the series of equation (14). Therefore, this coef- 
ficient is probably very small. Such a small mean latitudinal dependence of 
can readily be explained by the auroral heating according to equation (12) which 
compensates the respective (2,0) -term of the XUV heat input. Thus, we can 
estimate the ratio between the total heat inputs of corpuscular and XUV 
heating as 


^corp ^ 5 

Qyyv ^ 8 X 3.66 


0.17 


(15) 


A similar compensation occurs for the diurnal (l.l)-coefficients. However, 
here the effect is much smaller and yields a reduction of the total (l,l)-com- 
ponent by not more than 7% due to auroral heating. 


\ 








Secondly, the XUV heat input and the auroral heat input contribute to the (2.0)- 
component of the semiannual variation in the form 


I ^0 Qxuv 


- 0,2 X 3,66 


P“ cos 2 

Z a 


t > 


0.025 PO cos 2 


if we adopt the estimate of (15). This variation is very small. The observations 
within thermospheric heights indeed seem to support a rather small latitudinal 
dependence of the semiannual variation. The main contribution of the heat 
sources to the semiannual variations therefore appears to come from the (0,0) 
components of the auroral heat input (equation (12)) and from the heating due to 
wave dissipation (equation (13)). The sum of the three heat sources is then for 
the (0,0) component 

I 

<Qxuv ^ ^corp + Qdis>S = !• 17 Q^ov a - 0.029 cos 2 n t) 

+ (1 - 0.05 cos 2 t). 

With > Q°yy and < 0.17 , we conclude that the predominant 

generator of the semiannual variation is the energy due to wave dissipation al- 
though a small but significant fraction is also generated by XUV and auroral 
heating. 

According to an international convention (Hines, 1970) we shall call tidal waves 
all waves with periods of one lunar or solar day or a fraction of a day. These 
waves are obviously generated by those energy modes in equations (9) and (12) 
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with wave domain numbers m > 1. Planetary waves are wave modes with 
periods exceeding one day. They are therefore generated by the energy 
modes in equations (9), (12) and (13) with a wave domain number m = 0. Tidal 
waves will be discussed in the second part, planetary waves will be treated in 
the third part of this paper. 

4. APPROXIMATE GENERAL SOLUTIONS OF WAVE PROPAGATION 
In this section we want to solve simultaneously the equations of conservation of 
mass, momentum and energy in terms of eigenfunctions of the atmospheric dy- 
namic system taking into account the approximations and assumptions dis- 
cussed in section 2. In the following parts of this paper we shall use these 
solutions to discuss several tidal and planetary waves. 

According to our assumption in section 2 the quiet steady state structure of the 
atmosphere is related to the energy input of the (0,0)-components in equations 
(9), (12) and (13); 

I 

. ^*^WV ^oorp ^0 ® 

, I 

We are interested in any departure in time and space from that quiet atmospheric 
structure represented by the height- dependent mean values of temperature TJ, , 
density Pq * and pressure p^. The pertui-bated physical values may be written 
as temperature T, density p, pressure p, northerly wind u, westerly wind v, and 
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vertical wind w. The equations of conservation of mass, momentum and energy, 
and the equation of state are then presented for these perturbated amplitudes 
in spherical co-ordinates (r, B. \) 


— + di V V. 

Po 3 t hor 


3 w 
3 z 


w 

hT 


= 0 


+ 0.6 ^'u-2ncv-J_ 

3 t Po 



+ 


J_ 
^0 " 



3 V - „ 1 B / 

+ vv + 2 izcu- — — 77 
3 t ^0 ^ ^ \ 



_i_ lE = o 

r s 3 \ 


(16) 


g P + 1^ = 0 
3 z 


3 T 

^0 3 t ^ PQ 




L lA 

Po 




Here, for abbreviation, 


s = sin 6 *: C = cos 6 \ v^^^^ = (11, v). 


The factor 0.6 v in the second equation of (16) takes account of the latitudinal 
mean of the ion-neutral collision number u . 


26 


The heat source Q in (16) is the sum of solar XUV corpuscular heat input and 
dissipated energy and includes also a heat sink due to infrared cooling of the 
thermosphere. That cooling may be about 10% in magnitude of the XUV-heat 
input (Harris and Priester, 1962). But in view of the uncertainties in the adopted 
heating rate we shall not evaluate this heat sink but rather consider it as in- 
cluded in the effective XUV heat source. 


We develop the wave amplitudes of pressure p, temperature T, density p , and 
vertical wind w into series of the eigenfunctions of the system (16): 



(z)) 



il . m 1 a « t 

V 

(T^I'Wz)) 


(w^-f (z)) 


> ^ (0, z) exp [j (m \ + 2 f U t + s t)] ( 17 ) 


With 




(18) 


Furthermore, we define 


di v • V 


2 j f n H, 


hor 


V — I p"’ ' ^ _ 

■/ — ^ exp [j (m \ + 2 f fi t + s n t)] , (19) 

i—J hm, f " ® 

n 
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where we introduced the following abbreviations: 
f s n^)/2fi 

Cq = \'y gfip is the velocity of sound 

t 

f = w/ 2 n the Coriolis parameter 

(Ai is the angular frequency of the wave mode 
n is the angular frequency of the sidereal day 
is the angular frequency of the period of one year 
f is the equivalent depth of the mode of wave domain numbers (n,m,f) 
Hq is the scale height of the isothermal atmosphere 
6^^’' f (0, z) is the eigenfunction determining the latitudinal structure of the 
wave mode. 

For m ^ 1 and in the case that the dissipation terms can be neglected in (16), 



the eigenfunction 0'^'^ degenerates to the well Itnown Hough- function of tidal 
theory (Chapman and Lindzen, 1970). Contrary to this situation, our function 
f depends not only on co- latitude 6 but also on altitude z. Likewise, the 
equivalent depth hJJ' ^ which is real (either positive or negative) and which is constant 
within the lower atmosphere becomes height-dependent and generally complex 
within the dissipative atmosphere. This can readily bo seen from equation (16) 
if we neglect the Coriolis force in the second and third equation of (16) and in- 



troduce the dissipation factors from equation (3) and from equation 
(4). Then, the two equations of horizontal momentum in (16) become 
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1 




2 j Pq r 6 6 


V - 


1 Sp 

2 j Q Pq s r 9 \ 


with 


f, . f - i (0.6 Z„, . (2,.,)„) 


f. *(Zvi.U 


Taking for convenience 


(20) 


f ~ f ~ f, 

u V k 


(21) 


we obtain from equations (20) and (21) 


div V, 


hor 


2 j fk n r= 


A' p 


( 22 ) 


where A' is the horizontal component of the Laplace operator in spherical co- 
ordinates, The Laplace equation is satisfied by the spherical harmonic surface 
functitms 


YJJ = (0) e-’"^ 


(23) 


and yields 


A' Y" = - n (n + 1) Y" 


(24) 


Thus, the spherical functions become the eigenfunctions of system (16), and it 
follows from (18) and (19) 


Consequently, 


and 


f ^ pm 
n n 


Sm.f, 
n . n 



for n' = n 
for n' n 


4 f;^ f r2 

n (n + 1) g 


(25) 


(26) 


At a fixed height the eigenfunctions of (17) are a complete orthogonal system 
(Chapman and Lindzen, 1970). Therefore it is possible to develop any external 
heat source Q into a series of these eigenfunctions; 


Q.-:. 




J-n.f z) 

n iTi s f 


exp [j (m \ + 2 T n t + s n t)] 


(27) 


with 


.-- 1 ^ 


y 

X 0’"'^ 
TtTli f _ * ” 


In (27) we allowed for any angular frequency w of the energy source, not 
necessarily commensurable with the frequency of one sidereal day That 
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generalization may be important for the lunar tides and especially for the study 
of individual geomagnetic storm effects which contain a continual spectrum of 
periods (Volland and Mayr, 1971b). Moreover, it may be relevant for the study 
of tidal motions on the planets Venus and Mercui'y where solar day and sidereal 
day are quite diffei’ent. 

From equation (27) follows that 


B \ 


j ^ = j + s j 


(28) 


where we generally neglect the time derivative of exp [j s t] . 

Apparently, because of the height dependence of the eigenfunctions f , these 
eigenfunctions are coupled with each other in a realistic atmosphere. In this 
approximate treatment we shall neglect that coupling in order to obtain tractable 
solutions. That is, we assume 


1 






pm, f 

d z 


>'> 

1 

'd O’”’ f 

n 

1 

2l 

n 


n 

d z 

n 

3 z 

J 





In paper I we shall show that this assumption is reasonable at least outside the 
height region between 100 and 200 Km. 
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We are now prepared to derive from system (16) the equations for the height 
dependence of the coefficients p, w, and T, Taking account of the various as- 
sumptions and approximations in sections 2 and 4, furthermore using the ap- 
proximation equation (21i (which we shall abandon in the numerical study of 
paper I) we find from equation (16) the following system of ordinary differential 
equations of first order for p and w: 


with 


1 d e 


K e + h 



\ / 



2 ko H, n 


(30) 
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c.' -- 211. 

“ 2 c„ 


- f - j (Z|,^ from Equ. (2)) 
while the temperature amplitude is 


Tm. f - ^0 (l\ kK 

" 2 c„ 


_ i J"' f 


(31) 


Applying a well Imown matrix calculus (e.g. Volland, 1968), we find a phase 
integral solution of (30) (which is approximately valid as long as the elements 
of K in equation (30) change only slightly with height, and which is correct for 
constant elements of K^) to be 


c - J 
where 




f* ^ 



* / 

ko exp 

''o 

j 

I<o N d 

P“^ h d T exp 

j* 

N d 

^0 


c (Zq)( 32 ) 




\ expjj ( - j kg q + 1/2 Hq) d ^ 
t>o expjj (j q t 1/2 H^) d 




(33) 


are the upgoing (a) and downgoing (b) characteristic waves which are connected 
with the physical wave parameters of equation (30) by the transformation 


e = P c. 


(34) 
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The transformation matrix P is 



and P'^ its reciprocal. 


( 35 ) 


Moreover, the eigenvalues q can be found from the normalization of the co- 
efficient matrix K: 


(36) 


(37) 

The up- and downgoing waves are so defined that their amplitudes remain 
bounded in the direction of propagation. Thus, the imaginary part of the eigen- 
value of the upgoing waves is negative, and positive for the downgoing waves. 
Moreover, we laiow from the theory of gravity waves (Hines, 1960) that the 

‘Obviously, the subscripts and superscripts (™'^ should be added to the matrices e, c, h^ P, K, 

N as well as to the parameters q, F^, and We omit these indices for convenience as long as 
no confusion is to be expected. 


^-q-jA 0 


N - P’l KP : 


as 


q - j A> 


^ - A(a'"'f f j ,6"’ r. _ A 

n ' n ' n ' 


- 1 


>C' 'h 


vertical phase velocity of upgoing gravity waves is negative. Thus, the real part 
of the eigenvalue of upgoing waves is negative, too, and the terms a and /3 in 
equation (37) are both positive definite. That definition includes the two special 
cases of wave propagation within the nondissipative atmosphere in which 

There, it is p = 0 for the propagation modes, and a = 0 for the 
evanescent or trapped modes. Trapped modes with negative equivalent depths 
have an imaginary eigenvalue of /5 > 1. These waves can not transport wave 
energy in the vertical direction. Their wave energy leaks away vertically to 
such extent that the wave amplitude of the upgoing wave decreases according to 

a CO exp {- kg (/? - 1) A (z - z^)} (38) 

Thus, the adiabatic increase of the wave amplitude due to the exponential factor 
1/2 Ho in equation (33) is overcompensated by the leakage effect. We shall see 
in the following that generally within the dissipative thermosphere all wave 
modes have imaginary terms ,6 > 1, but real terms a different from zero. 
Therefore, we shall call waves of that type quasi- evanescent modes. They ful- 
fill the radiation condition as well as the condition that the kinetic energy density 
is bounded at z Our eigenvalue equation (37) of course degenerates as it 
should to the eigenvalue of the nondissipative atmosphere for = f|^ - f (Chap- 
man an.^. Lindzen, 1970). 
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5. SPECIAL SOLUTIONS OF THE GENEIiATION AND PROPAGATION 


OF ATMOSPHERIC WAVES 

We now want to discuss two special solutions of equation (30) which have some 
relevance to the generation and propagation of waves within the atmosphere. 

We shall restrict ourselves in this section to an isothermal atmosphere with 

constant dissipation factors Z; that is, constant elements of the coefficient 

matrix I£ in equation (30). A more quantitative treatment including realistic 

temperature profiles and vertical profiles of the dissipation factors will be 

given in paper I. 

a. Infinitely extended isothermal atmos phere 
We first discuss the simplest case of wave propagation within an isothermal 
atmosphere infinitely extended in the vertical direction which shall simulate wave 
generation within the thermosphere. We assume an external heat input with the 
vertical profile (see equation (5)) 


jm, 

n 




)/2H„ 


fo r 


^ z ^ Zj 


(39) 


Outside that region the heat input is assumed zero. For the boundary conditions 
we adopt the simple radiation conditions namely that waves generated within the 
region z ^ < z < z ^ by the heat source of equation (30) can only leave that region 
through Zq and Zj, respectively. This implies that 


« (Zq) = b (Zq) = 0 
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The explicit solution of (30) then becomes 


a (z) ~ 0 


b (z) - b (Zq) exp [j (q ~ j A) (z - z^)] 


a (z) = (1 - exp [- j (q + j e A) (z - z^)]) 


b (z) = {1 - exp [j k (q - j e A) (z - Zj)]} 


fo r z z„ 


(41a) 




a (z) = a (Zj) exp [- j k^, (q + j A) (z ~ Zj)] 
b ( z) =0 


for Zj < z 


(41c) 


with 


j J (f H A) 

(K - F^) (q + j e A) 


j J (f F^ + j A) 
fh fFb - FJ (q - j 6 A) 


The physical wave parameters are according to Equations (31) and (34) 


W 


m, f 


"•0 

> m, f 


• = a + b 


■= F^ a + b 


( 42 ) 
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T 2 f 

*0 ^ ‘h 


[2 K ( f + j A) a T 2 /< ( f F^ + j A) b - j J] 
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where 


.< Q*"' f 

T „ xm, f _ 


from equations (39) and (27). 

With the boundary conditions (40) we implicitly neglected the reflection of waves 
(generated within Zq < z < Zj) at the earth's surface. In an exact treatment one 
would have to add to the upgoing wave in (41) the wave 

a (z) - b (Zq) exp {- 2 j q » j k^^ (q .. j A) (z - z^)} 

which is the downgoing wave of equation (41a) reflected at the earth’s surface. 
For the trapped modes this reflected wave is insignificant at z ^ 100 km. For 

the propagation waves, a is only significant in the immediate vicinity of the 
lower boundary z ^ and it is entirely negligible when compared with the wave 
generated below z^ (see section 5b). Thus, our solution (41) describes suf- 
ficiently well the generation and propagation of waves within the thermosphere. 


From the approximated equivalent depth at thermospheric altitudes in equation 
(26) and from the equation of the eigenvalue q in (37) we estimate that a wave 
mode becomes quasi-evanescent in a region where 


38 


which is certainly valid above about 200 km altitude. There the imaginary part 
of the eigenvalue in (37) becomes B > 1, and we find from equation (41b) that 
the waves approach the asymptotic values 

a - G 1 

A I 

I for 00 

b^G, 

if the heat source equation (39) is extended into infinite. Then, the wave ampli- 
tudes of equation (42) approach the asymptotic values 


f J (2 Hg/h ^ 1 - 6 ) ^ J [n (n + 1) - 2 (1 T e) f 

A[(62„1) f,^x4^.fH^/h] 2A[(..2_1) f^ ^ + k n (n + l)j 


pm, f 
^ n 

•pm, f 


Fa Ga + F, G, = 


j J (1 ^ 


j J (1 + f, S'2 y 


(45) 


(e2 _ 1) + 4^ f H/h [(.2 _ 1) ^2 n , 1)] 
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and from the last Equation (16), 
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2 Po 
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^ - 
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Here, the last equations on the right hand side of equations (45) have been cal- 
culated with help of the approximate equivalent depth of equation (26). 


From equations (20) and (26) we find the assymptotic wave amplitudes of the 
horizontal winds 


V 

c~ 


m. f d P"’ 

"" “ {j(m\+ 2 f fl t + s t)} 


111) . f pi’i 

" j Q, ~T"^ exp (j(m\ + 2 ? n t + s O t)} 


(46) 


” s i n 0 
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J (1 + o 


2 [(e2 _ 1) y ^ ^ n (n + 1)] 


Since f,, and f ^ tend to become constant \vith altitude (see section 2) the relative 
amplitudes of the waves also approach constant values if e= 1, and we obtain 
for z -• 00 


P/Po 


^■p/po^ 


2U 

X n (n + 1) 


U' C„ 00 V''C^ 00 


n (n 4 1) 


T/Tq - 0 

w/Wq (1 _ 4 f ^2 y ^/n (n + 1)) 


(47) 
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That means, the thermosphere behaves like a low pass filter which suppresses 
the pressure, density and horizontal wind amplitudes of the wave modes pro- 
portional to l/n(n + 1), The same is the case for the vertical wind for low wave 
domain numbers n, while the vertical winds tends to become independent of n for 
large numbers n. The temperature amplitude approaches zero at great heights 
which is not consistent with the observations. However, it will be seen in the 
following parts of this paper that this decrease with altitude occurs very slowly and 
that the assymptotic value is only reached at heights of several thousands of 
kilometers. Below 400 km, the temperature remains nearly constant. As al- 
ready mentioned before the temperature amplitude is significantly affected by 
the neglect of heat conduction waves. In paper I it will be shown that due to 
heat conduction waves the temperature amplitude slightly increases toward a 
finite assymptotic value at high altitudes. 



We can simulate this influence of heat conduction waves in our simplified model 
by introducing in equation (5) a value of 


I e ' 0.75 

1 


(48) 


which corresponds to a heat input that decreases with altitude somewhat slower 
than the mean pressure. Then we find the ratio between temperature and density 
from (45) 
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T/To 

P/P^ 


1 •- € 

r- i — : '-0.14 
1 - e 


which is slightly smaller than in the Jacchia- model (Jacchia, 1964). 


Taking the value of e from equation (48) and assuming plausible physical param- 
eters for the thermosphere we find from equation (46) a formula for the effective 
assymptotic magnitudes of pressure, density, temperature and horizontal winds 



Qm. f 
^ n 

0.4 Y n (n + 1) 


(49) 


This formula is valid within the height range between about 350 and 400 km and 
is a very convenient half empirical formula (Volland and Mayr, 1971b) for tidal 
and planetary waves at thermospheric heights as verified by the theory (see 
paper I and the results of the following parts of this paper). 

b. Wave generation within the lower atmosphere 

We now want to consider wave generation within the lower atmosphere which 

is important for the study of the tidal propagation modes. The relative amplitudes 
of tidal propagation modes with zero attenuation (/? = 0) increase with height 
according to exp (Az) within the nondissipative atmosphere. We therefore expect 
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significant wave amplitudes for the upgoing waves of those modes at the base of 
the thermosphere which we have to add to the upgoing waves generated within 
the thermosphere. The lower boundary condition at the earth's surface is that 
vertical winds should disappear. Then, taking a heat source of the form of 
equation (39) within the height range between 0 ^ z ^ z ^ we obtain from 
equation (32) 


n ( Z) •' 1 - rxp - j k„ (q j • A) 


G|, 1 - PXP I - j (q - j , A) z^; - csp I- j (q . j A) z. 


for 0 < z ? Zq 


(50a) 


b (z) r. G,, (1 „ exp [j (q - j e A) (z - z,,)]} 
and 


a (z) r: a (Zq) exp [- j (q ^ j A) (z - z^)] 

for z < z 

b(z) = 0 


(50b) 


Because of ^ 0, it is at the earth's surface (z = 0) 


a = - b 


for z = 0 


b ~ G, 

b 


Thus 


(51) 





j(fF3 + jA) 
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Pq f (q ■' j e A) 


for z = 0; f,, = f 
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At the base of the thermosphere (z = Zg ) it is for propagation modes: 

a (Zq) '' - j ><o *■ j ^ ^’^P t- j *<o + j ^0^'^ 

for5<l, (52a) 

b (Zg) .. 0 

and for evanescent modes: 


a (Zq) 

for /S > 1. (52b) 

b (Zg) - 0 

The observed atmospheric wave modes are the sum of the wave modes generated 
by the energy input within the lower atmosphere (equation (50)) and the energy 
input within the thermosphere (equation (41)). 

6. CONCLUSION 

In this paper we developed a three-dimensional model of thermospheric dynamics 
in terms of the eigenfunctions of the atmospheric system. These eigenfunctions 
or wave modes are excited by solar heating from XUV-radiation from particle 
precipitai.'ou and joule heating in the auroral zone during geomagnetic disturbances 
and from energy dissipation of waves from the lower atmosphere. We determined 
formulae for both solar heat sources in terms of the eigenfunctions of the 
atmospheric system. That series contains tidal components depending on local 
time and planetary components depending on seasonal time. We estimated the 
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relative importance of auroral heating when compared with solar XUV-heating 
and found that the average auroral heating contributes not more than about 15% 
to the total solar heat input. 

Approximate analytic solutions for the generation and propagation of atmospheric 
waves within the dissipative thermosphere excited by the solar heat input have 
been derived. It was shown that the eigenfimctions which are the Hough- functions 
within the non-dissipative lower atmosphere, change into the spherical surface 
functions within the dissipative thermosphere. The amplitudes of density and 
horizontal winds for the various wave modes are shown to decrease proportional 
to 1/n^ where n is the zonal wave domain number of the spherical harmonics. 
Therefore, only wave modes with low wave domain number n are significant at 
thermospheric heights. 

In two further parts of this paper we shall discuss in detail the characteristics 
of various tidal and planetary wave modes, especially the height dependence of 
the eigenfunctions and equivalent depths and the change in the latitudinal structure 
of the wave modes. In an additional paper (Volland and Mayr, 1971a) numerical 
full wave calculations are carried out to develop a more sophisticated and 
quantitative thermosphere model. 
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